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PRESENTAZIONE: Si considera un sistema bidimensionale della dinamica dei gas introdotto nel 
1970 da Ya. Zeldovich per descrivere la formazione della struttura di grande scala dell' universo. II 
sistema si rivela come qualcosa di intermedio tra un sistema di equazioni differenziali ordinarie e un 
sistema iperbolico di equazioni alle derivate parziali. La caratteristica principale e' la nascita di singo- 
larita': discontinuita' della velocita' e funzioni delta di vario tipo per la densita'. Si da' una descrizione 
rigorosa delle soluzioni generalzzate in termini di misure di Radon e si ottiene una generalizzazione 
delle condizione di Rankine- Hugoniot. Sulla base di tali condizioni si mostra che la rappresentazione 
variazionale delle soluzioni generalizzate, valida nel caso unidimensionale, non vale in generate nel caso 
bidimensionale. Si ottiene anche un sistema unidimensionale non banale non strettamente iperbolico per 
la descrizione dell' evoluzione delle all' interno dell' urto. 

ABSTRACT: The paper deals with the 2-D system of gas dynamics without pressure which was 
introduced in 1970 by Ya. Zeldovich to describe the formation of large-scale structure of the Universe. 
Such system occurs to be an intermediate object between the systems of ordinary differential equations 
and hyperbolic systems of PDE. The main its feature is the arising of singularities: discontinuities for 
velocity and 5 -functions of various types for density. The rigorous notion of generalized solutions in terms 
of Radon measures is introduced and the generalization of Rankine- Hugoniot conditions is obtained. On 
the basis of such conditions it is shown that the variational representation for the generalized solutions, 
which is valid for 1-D case, in 2-D case generally speaking does not take place. A nontrivial 1-D system 
of nonstrictly hyperbolic type is also obtained to describe the evolution inside the shock. 
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1 Introduction. Basic definitions. 



This paper studies the shock waves for the system of 2-D pressureless gas 
dynamics 
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where (u, v) has the physical meaning of velocity vector, g has the physical 
meaning of density. The system (HD is nonstrictly hyperbolic system of con- 
servation laws which has three coinciding characteristic fields and incomplete 
set of eigenvectors. Due to these properties the shock waves develop strong 
singularities in the density which are of type of ^-functions on the surface. 
So ([l]) is relevant to describe some process of concentration of matter. We 
will consider the Cauchy problem for (p]) 



Q(0,x,y) = go(x,y) > 
u(0,x,y) = u (x,y) 
v(Q,x,y) = v (x,y) , 



(2) 



where go, uq, Vq are piecewise C 1 (1R 2 ) functions and will be taken in the special 
form (see below) to produce locally the single shock front. 

For the smooth functions the system ([[]) is equivalent to the following 
system of equations 



dg d(gu) d(gv) 
dt dx dy 

du du du 

dt dx dy 







(3) 
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The last two equations in the system (H) constitute the so-called inviscid 
2-D Burgers equation, which was proposed (but in 3-D case) by Ya. Zeldovich 



24 1 to describe the formation of the large scale structure of the Universe. 



Further this approach, in particular the consideration of the whole system 
(0), was developed in the consequent papers (see, for example, [[HJ, fl2|| , p2 



and the references therein) from the physical point of view. But if we 
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are dealing with the laws of conservation of mass and momentum it would be 
more convenient from the mathematical and physical also points of view to 
investigate the system in divergent form (|I|) especially in the case when one 
has the developing of such singularities as shock waves. For 1-D variant of the 
system flU) the particular generalized solutions were constructed in |2| (there 
also were proposed some numerical schemes), and the existence theorem in 
the sense of Radon measures for wide class of initial data was proved in |§ , 
TU| , [|11|] , see also [[J , |4j] for the uniqueness results in the framework of duality 



solutions. In ||, [|Kj the generalized solutions were constructed by the aid 
of variational principle for appropriate Hamilton- Jacobi equation (see |T3 



fll4f , fl5fl , [|l7j1 for the variational principle to one quasilinear equation of the 
first order). 

There were number of attempts to construct the generalized solutions to 
system (fj) by analogy with 1-D case and with the aid of variational principle, 
which can be written for 2-D inviscid Burgers equation and then general- 
ized through corresponding Hamilton- Jacobi equation. Further in this paper 
we will show that in general this is impossible. The problems for multi-D 
which involve different types of singularities can expose rather complicated 



behavior, see, for example, |TBL M, p5f. Here we investigate the 'simplest' 



degenerate case of 2-D system when the characteristics can be calculated ex- 
plicitly. Nevertheless the generalized solutions to system (|I|) (see Definition 
1 below) exhibit the different behavior compared to ones for the single first 
order 2-D equation of Hamilton- Jacobi type. The main ideas of the present 



paper were outlined in El, |20] 



As it has been already said for smooth functions the system (p]) is equiv- 
alent to (^) and for small enough values of time and smooth initial data the 
characteristics method can be applied. So the solution to the problem (P, 
(0) reads 

Q[t,x,y) = eo(a,b)— 

d(x,y) (4) 

u(t, x, y) = u(0, a, b) ; v(t, x, y) — v(0, a, b) , 

where functions a(t, x, y), b(t, x, y) are determined by the equations 

x = a + tu(0, a, b) ; y — b + tv(0, a, b) . (5) 

But it is well known that the characteristics of the system ([3]) can intersect 
themselves for finite time even for infinitely smooth initial data. Hence one 
finds the formation of singularities: discontinuities for the velocity (u, v) and 
5-functions for the density g which correspond to concentration process in 
some points or along some curves. So it is necessary to introduce the notion 
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of generalized solution to the problem ([I]), (^]) which is natural to formulate 
in terms of Radon measures because of the type of singularities. 



DEFINITION 1 Suppose (P t (dx,dy), I t (dx,dy), J t (dx,dy)) are the families 
of Radon measures defined on Borel subsets of M? , weakly continuous with 
respect to t and such that P t > 0, and I t , J t are absolutely continuous with 
respect to P t for almost every fixed t > 0. Let us define the vector function 
(u(t,x,y),v(t,x,y)) as Radon-Nykodim derivatives 

Then (P t , I t , Jt)) will be called the generalized solution of the problem (fj^ (^), 
iff: 1) for an arbitrary functions f,g,h G Cg(lR 2 ) (the space of continuously 
differentiable functions with compact support) and < t\ < £ 2 < +oo 



f(x,y)Pt 2 (dx,dy) - J J f(x,y)P tl (dx,dy) = 

df 

—(x,y)I T (dx,dy)+ (7) 



g{x,y)It 2 {dx,dy) - J J g(x,y)I h (dx,dy) = 
dg 

—{x,y)u(T,x,y)I T (dx,dy)+ (8) 

h(x,y)J t2 (dx,dy) - J j h(x,y)J tl (dx, dy) = 
dh 

— {x, y)u{r, x, y)J T {dx, dy) + (9) 



where J J stands for the integration with respect to whole 
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2) in the weak sense as t — > +0 

P t — > f?(0, a, b)dadb ; J t — > £>(0, a, 6)w(0, a, b)dadb ; 
J* ~~ * £?(0, a, 6)f (0, a, b)dadb . 

Concerning the system (|3|) there is well-known Hopf-Cole's representation 
13 1 for the solution with potential initial velocity. This representation can 



be derived with the aid of infinitesimal viscosity method (i.e. the adding of 
the second order operators with small parameter e to the right-hand side of 
the last two equations of (|3|) and then letting e tend to zero) and allows to 
determine the location of singularities. Namely 

(u,v) = V( XlV) V(t,x,y) , (10) 

where 

. r ry , (x - a) 2 (y - b) 2 

a, 6 Zt It 

So is the potential of initial velocity. The representation (0) is valid in do- 
mains of smoothness of \1/ and gives the location of the set where singularities 
arise. That is the singular points of \I/ will be such points (t, x, y) that the 
global minimum of the expression in the braces with respect to (a, b) is at- 
tained more than in one point. For 1-D case the representation (pi]) in case 
of constant initial density gives us the way to construct the generalized solu- 
tion to the 1-D analogue of the problem ([]]), (|2|), see 0, [K| (in these papers 
one can also find the generalization of (|10"D to nonconstant initial density). 
Unfortunately in 2-D case the generalized solution in the sense of Definition 
1 to the problem (|l|), (0) for constant initial density can not be constructed 
via(0). 

Our study is organized as follows. In § 2 we formulate the conditions on 
the initial data (0) which allow to conjecture the development of shock and 
then assuming that such shock exists and satisfies the problem ([]]), (fj) in 
the sense of Definition 1 give the formulas for the velocity along the shock. 
These formulas are the generalization of Rankine-Hugoniot conditions for 
hyperbolic conservation laws. In § 3 it is shown that in the case of potential 
initial velocity these formulas determine the singularity surface which in gen- 
eral does not coincide with one determined by the representation fllCf) . As 
it is shown in § 4 one can obtain certain 1-D quasilinear hyperbolic system 
of PDEs for the motion inside the shock. But this system does not sat- 
isfy Friedrichs' symmetrizability condition, so it creates difficulties to prove 
the global existence theorem. An interesting heuristic procedure to guess 
the system is deferred to Appendix. Finally in § 5 one investigates some 
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simplest solutions for derived system in case when g , uq, vq are piecewise 
constant but the data along the initial shock front are rather arbitrary. Even 
with such trivial external field one encounters with complicate behavior of 
the generalized solutions. It is shown that there are nontrivial cases when 
the solution can be defined explicitly and other cases when the correspond- 
ing Cauchy problem is reduced to equation of type P tt = const ■ P x . Such 
equations are ill-posed in the spaces of functions of finite smoothness. 

Further the letter subscripts (except Y, , f and '£') will denote appropriate 
derivatives. 

2 The propagation of the shock front. 

Let us denote the coordinates in R 2 x {t = 0} as (a, b). Suppose there is the 
C 1 curve Q = (A(l), B(l)) in the plane (a, b), I is a parameter, such that, for 
definiteness, Ai > 0, Bi > ; Af + Bf ^ . Suppose that Q can also be 
found as the solution of the equation G(a, b) = 0, G G C 1 (R 2 ). Consider the 
following initial data (|2|) 

g (a, b) = £u + (g + - £>_) H(G(a, b)) 

u (a,b) = m_ + (u+ -u-) H(G(a,b)) (11) 

t>o(a, b) = V- + (v + — V-) H(G(a, b)) , 

where H is the Heaviside function, i.e. H(6) = for 9 < 0, H (6) = 1 for 9 > 0; 
and the following conditions hold 

I) u_(o,6),i;_(a,6) G C l (<F) ; u+{a,b),v+{a,b) G C l {G+) , 
here Q~ = {{a, b) : G(a, b) < 0}, G + = {{a, b) : G{a,b) > 0}. 

II) Aiv_ - Biu_ < and A t v + - B t u + > on Q. 

REMARK 2.1 Let us note that the condition I) reads that for every point 
(a*,b*) = (A(l*),B(l*)) G Q there exists some domain Q C M 2 x [0,T(Z*)] ; 
T(l*) > 0, (a*, b*) G Q f]{t = 0} where the transformations (j = +, —) 

x = a + tUj(a, b) 
y = b + tvj (a, b) 

are nondegenerate for every < t < T(l*), the characteristic lines (W with 
u(0,a,b) = v(0,a,b) = V- never cross themselves, and the characteristic 
lines with u(Q, a, b) = u + , v(0, a, b) = v + never cross themselves. 

But the condition II) reads that for < t < T(l*) there will be another 
domain Q 1 C Q where exactly two characteristic lines issued from different 
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sides of Q will come to the same point (among these points the shock surface 
will form itself). So these conditions are natural for the propagating shock 
front to exist locally up to time moment T(l*) > 0. 

The following definition is essential for our analysis to elicit shock fronts 
with 'good' behavior. Let us denote through Qi the maximal Qi-domain 
from Remark 2.1 which is valid for 1 = 1*. 

DEFINITION 2 The shock front S will be called stable iff 

S P|{t = r} C Qi f]{t = r} as < r < T(l*) . 

It is natural to seek the solution to the problem ([I]), (JTTJ) in the form (see 
Fig. |I] below, Eulerian representation) 

g(t, x, y) = g- + (f?+ - Q-) H(S) + P t 5(S) 
guit, x, y) = Q-U- + (g+u+ — Q-U-) H(S) + ItS(S) (12) 
gv(t, x, y) = g_v_ + (g+v + - Q-vJ) H(S) + J t S(S) , 

where g j (t,x,y),u j (t,x,y),v j (t,x,y),S(t,x,y) E C 1 (R + x R 2 ), (j = +,-) 
and satisfy ([I]) in classical sense; S(t, x, y) = represents some surface S 
in R + x R 2 and S f] ({0} x R 2 ) = Q; P t ,I t , X e C^S); H is the Heaviside 
function mentioned above, 8 is standard Dirac measure. 

If the solution to the problem (|]), (|ll|) exists in the sense of Definition 
1 then by formulas ([|) there is defined the velocity vector U = (u, v) on 
the surface S and from each point of Q one can draw the integral curve 
of U on S. Suppose (x s (t,l),y s (t,l)) is the corresponding parametrization 
of S (I is a parameter along Q) and d(x s ) jdt = u, d(y s ) jdt = v. Then 
according to formulas ([D and conditions I), II) one can define the functions 
aj(t,l),bj(t,l) G C 1 (iS), (j = +, — ) which correspond to the initial positions 
of characteristics which come to the same points of the surface S from left 
and right. 

Further let us define the map (see Fig. [I] below) 

C t : (a, b) -> (x,y) (13) 
in the following way. Let us issue the characteristic line (keeping in mind 



(|lTl) ) from the point (a, b) 

x(t) = a + Tu (a,b) ; y(r) = b + rv (a,b) (14) 
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and consider the time To of the intersection of line ( fn|) and S. Then let us 
take C t (a, b) = (x(t),y(t)) if r > t and £f(a, 6) = l ),y s (t, l )) if r < t, 
where l is defined from the condition (a;(r ), y(r )) = (x s (r , l ), y s (T ,l )). 

Now we are able to formulate the generalization of Rankine-Hugoniot 
conditions. 

THEOREM 2.1 The generalized solution in the sense of Definition 1 to the 
problem (fjj, ( fil| j in the form exists iff the following formulas are true 

u = I t /P t ; v = J t /P t (15) 



o 

t 



[Qo ((a+)r(b+)i - (6+) T (o+)i) - g ((a_) T (6_), - (6_) T (a_)j)] dr 

h= I [Qo Uq (( a +)r(6+)« - ( & +)r(a+)«) - £ U (( fl -)r(&-)« ~ (6_) T (a_)i)] ^ 

(16) 

Jt = I [qqVq ((a + ) T (6+)j - (6+) T (a + )j) - £ v ((a_) T (6_)j - (6_) T (a_)j)] dr, 



o 



w/jere ^ = ^(a,-, &,•)> w£ = m (%, &,•), v 3 = v (aj, bj), (j = +,-). 

PROOF. Suppose there are some family of domains D(t) with oriented 
boundaries dD{r) in the plane (a, b), D(ti) C Dfo) for < t\ < r 2 < T. 
Suppose 0D{t) is closed curve (a(r, /), 6(r, /)), / is a parameter along the 
curve, a, 5 G C 1 ([r 1 ,r 2 ] x JR.). Then it is easy to check that the following 
formula is valid 

~j~ J J ydadb = J J ^-dadb+ j> ip(a T bi — b T a t )dl , (17) 

D(t) D(t) <9D(t) 

where <p G C 1 ([n, r 2 ] x U re[ri , T2 ] W)) ■ 

Denote through II some circle in the plane (a, b) such that II x {t} f] S t ^ 
0, S t = S f]{t} xR 2 . Denote through £L(t), £>+(*) (Consult with the Fig. §) 
the domains in the plain (x, y) to which at time t the characteristics come 
from (a_,6_), (a+, 6+) respectively (note that -D_(£) |J D + {t) = II x {t} 
and D_{t) f] D + {t) = S t ). Further let us consider the right hand side of 
integral identity (|7|). Take / G Cq(1R 2 ) such that supp f C II. Suppose 
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Eulerian Lagrangian 




a-(t,l),b-(t,l)) 
(a + (t,l),b + (t,l)) 

Figure 1: The map £ t for fixed t; I is a parameter along S. 

D*(t) = C^iDjit)), (j = +,-), D s (t) = C^iSt). Define the function 
ft(a,b) such that / t *(o, b) = f(£ t (a,b)). Then one has 



jdtlj^JJ (f x u + f y v)P t (dx,dy) + j (f x u + f y v)P t (dl) 

h \j=+- D u) St 



/*{ £ 1 1 



f x (a + tu , b + tv ) g (a, b)u (a, b) + 
f y (a + tu ,b + tv ) Q (a,b)v (a,b)^dadb + J f x (x s (t,l),y s (t,l)) u(t,l) 

df? 



+ 



St 



f y (x s (t,l),y s (t,l))v(t,l) P t (dl)\ = dtl 



ti 




dt 



■Qo(a, b)dadb + 



j^L( x ^t,l),y s (t,l))P t (dl)\ = fdt\Y.J t j j ft*Qo(a,b)dadb 

St ' ti U=+- D *( t ) 

J f(x s (t,l),y s (t,l)) 



s t 



Oo(a-,b-) 



d(a_)d(b_) d(b_)d(a_ 



y±( x °(t,l),y s (t,l))P t (dl)} , 



dt dl dt dl 
dl + 



s t 



using formula QT7|). 

Consider the right hand side of integral identity @, here g is playing the 
role of /, 



J dtlj2 J J (9xU + g y v) I t (dx,dy) + J (g x u + g y v) I t (dl) 
ti \i=+- Dj{ t) St 



*\ E 




g x (a + tu , b + tu ) £ (a, b)u 2 Q {a, b) + 



^ (a + tw , & + tv ) g (a, b)u (a, b)v (a, b) 



dadb + 



g x (x s (t, I), y s (t, /)) u(t, I) + g y (x s (t, l),y s (t, /)) v(t, I) I t {dl) 



t2 , 




9t 



Qo(a } b)u (a } b)dadb + / — ^ (x s (t,l),y s (t,l)) I t (dl) 



dt 



S, 



* E 




g^Qo(a, b)u (a, b)dadb — 



D*At) 



g(x s (t,l),y°(t,i)) 



g u 



<9(a_)<9(fr_) _ <9(fr_)d(a_) 



at a/ at a/ 



dt 



x s (t,l),y s (t,l))I t (dl) 



using formula (|TTD . 

Consider the right hand side of integral identity (|9]) , here h is playing the 
role of /, 



*2 



fl 



;'=+,- 




(/ia-u + /i y w ) J t (dx, dy) + / (/ij-u + hyv) J t (dl) 



DM) 



St 



t-2 



* E 




/i x (a + tu , 6 + tv ) g (a, b)u {a, b)v {a, b) + 



D'At) 
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hy (a + tuo, b + tv ) Qo(a, b)vl{a, b) 



dadb + 



/ 

Si 

t 2 



h x (x s (t, l),y s (t, I)) u{t, I) + h y {x s {t, l),y s (t, I)) v{t, I) J t (dl) 




dh* f dh* 1 

t -g (a,b)v (a,b)dadb+ —^(x s (t,l),y s (t,l))J t (dl)} 



dt 



dt 



D*At) 



J dt\ 7tt f f K&o( a ib)vo(a,b)dadb 



D*At) 



h(x s (t,l),y s (t,l)) 



s t 



0o v o 



d(g_)<9(6_) _ d(b-)d(a-) 
dt dT dt dT 



0o v o 



d(a + )d(b + ) d(b + )d(a + ) \ 
dt dl dt dl ) 

again using formula flTT|). 



dl + 



dh* 

-j£(x s (t,l),y s (t,l))J t (dl) 



□ 

THEOREM 2.2 For every I and < t < T(l) the following formulas are 
true 

t 

J [x s (t, I) - a+(r, I) - tu+(r, I)] g {a + , b + ) ((o+) T (6 + )j - (6 + ) T (o + )j) dr = 

[x s (t,l) -a-(r,l) -tUo(r,l)] g (a-,bJ) ((a_) T (6_)j - (6_) r (a_) { ) dr ; 

(18) 

[y s (t, I) - b+(r, I) - tv£(r, I)] g (a +1 b + ) ((a+) T (6 + )j - (&+) r (a+)i) dr = 

[y s {t, I) - 6_(r, /) - too (r, I)] g {a-, 6_) ((a_) T (6_), - (6_) T (a_),) rfr 
PROOF. Let us take formulas (p~5| ) and write them in the following way 

(X S ) T Pt=T = h=T ] (y S ) T Pt=T — Jt=T ) 
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t 







Now let us integrate these equalities with respect to r from to t and then 
integrating by parts in both sides one obtains 

t 

J (x s (t)-x s (r))x 
o 

[Qo ((a+) T (6+)j - (6+) T (a+)j) - g ((a_) T (6_)j - (6_) T (a_)/)] dr = 
(t - r) x 

[go M o (O+MM; - (6 + ) T (o+)j) - g u ((a_) T (6_), - (6_) T (a_),)] 

i 

(y-(t)-y'(r))x 

[ffo (( a +)r(&+)i - (6+) T (a+) { ) - £ ((a_) T (6_)j - (6_) T (a_)i)] rfr = 
(t-r)x 

[ffo u o" (( a +)r(Mz - (6+) T (o+)i) - Qo v o (( a -)r(b-)i - (6_) T (a_)j)] dr. 
Then applying formulas flT4]) we come to the assertion of Theorem 2.2. 



o 



o 



□ 

REMARK 2.2 Let us mention that expressions ftH\ ) look similar to the ex- 
pressions for the adhesion principle in JjQj, / fiZj/ . So all characteristics initially 



started at points (aj(r,l),bj(r,l)), j — +,—, r < t, I is fixed, will 'concen- 
trate' in one point at time t. 

3 Comparison with the variational represen- 
tation. 

In this section we assume that Qo(a, b) = 1. Now let us compare the formulas 
flIE|) for shocks with the restrictions on singularity surface which can be ob- 
tained through the variational representation (|T0|) in case of potential initial 
data. Suppose !F—(x v (t,l),y v (t,l)) is the singularity surface obtained from 
(flOl) . Then T can be represented as the solution of the equation 

F(a+,b+;t,x,y) - F(a-,b-;t,x,y) = , (19) 
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where 

F = S (a,&) + — ^— + — 5;— , 
and ctj(t, x, y), bj(t, x, y), (j = +, — ) are denned from the system of equations 
F a (a, b; t, x, y) = ; F b (a, b; t, x, y) = . 

THEOREM 3.1 Suppose the initial velocity vector in (|77[) is a potential vec- 
tor and the potential So (a, b) has the form 

S (a,b) = , asb> -ef(a,b) 
So(a,b) — b + ef(a,b) , as b < —ef(a,b) , 

where e > is sufficiently small and the function f G C°°(R 2 ), /(0,0) = 0, 
satisfies the conditions: there exists such neighborhood of the point (a = 0, b = 
0) that /aa(O) 0) ^ and the curve determined by the equation b+ef(a, b) = 
is monotone for every sufficiently small e > 0. 

Then the formulas j^TSj ) and / fJ^ j determine different surfaces. 



PROOF. It is easy to infer from that as e > is sufficiently small in 
the neighborhood of the point (t — 0, x\ — 0, x 2 = 0) the singularities surface 
has the form 

y = l-ef(x,- t -)+o(e). (21) 

Now let us find the equation of singularities surface on the basis of for- 
mulas fllSl) . As e = the solution of (|i~5|) looks as follows 

x(t, I) = a_(t, I) = a+(t, l) = l ; y(t, I) = 6_(t, /) = -6+(t, Z) = t/2 . 
As £ > let us seek the solution of the system (|15D in the form 

x(t, I) = Z + e x(t, Z) + o(e) ; Z) = t/2 + ey(t, I) + o(e) . (22) 
Then taking into account the relations 

x(t,l) = a-(t,l) = a + (t,l) +tef a (a + ,b + ) 

y(t, I) = b-(t, I) = b+(t, l)+t(l + ef b (a+, b + )) , 

one obtains P t , I t , Jt and finds 



dt 2t 



t t 

y(t, I) = -f(l, 0)-\J rfaadr + ±Jf b d T +±J r 2 f aa dr 



(23) 
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where the values of the derivatives of the function f(a,b) are taken at the 
point (/,— r/2). But functions (fffi) must satisfy (|2l| ) which in view of (|23|) 
is equivalent to the following identity 

t t 

f(l, -t/2) - f(l, $)-\J ( r " y) f~ dr + ~Jhdr = 0, 



where the values of the derivatives of the function /(a, b) are taken at the 
point (I,— t/2). Further, multiplying by t then three times differentiating 
with respect to t and substituting t = one obtains f aa (l,0) = 0. The last 
identity contradicts to the conditions on function f(a,b). 

□ 

Differentiating ( |l9l) with respect to t and changing notations from t to r 
one gets the following condition on the surface J- 

(a + - a_) (x*) T + (6+ - 6_) = 5 (a + , 6+) - 5 (a_, 6_) , (24) 



THEOREM 3.2 Suppose one has the potential initial data (|?~|). T/ien the 
surfaces which are defined by fiH\) and coincide iff the following relation 
is true 

t t t 

(a+ - a_) | J i*dr + (b + -b-)±J fdr = (S + -S.)± J p*dr , (25) 



where 



_ /, / , v + T / \ (6+)t + (&-) 

p ={b + - 6_) (a+ - a_) — 



(a + 


.)r + (a-)r 




2 




)r + {3-) T 




2 




f)r+(5-) T 


2 



f = (6 + - 6_) — - - (S + - 5-_) 
j* = (a+ - o_) (g+)r ± ( ^ )r - (S+ - 5_ 



2 

(Mr + (6-)r 

2 

(o + ) T + (a_) T 



PROOF. It is easy to see that because of our construction of the shock 
surface x s , y s , keeping in mind the assumption that x s = x v , y s = y v , one 
has 

X s = a_ + r-u_(a_, bJ) = a + + ru + (a + , b + ) 
y s = a_ + rt>_(a_, 6_) = a + + Ti>+(a + , 6 + ) . 
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and the relations ( p~5| ) at least locally determine the shock surface (if it exists). 
So (24) occurs to be an additional relation and it is consistent with (|15"D only 
if it follows from them. 

Hence the condition (|24T) can be rewritten as follows 



{a + -a-)I t + (b + -b-)J t = {S {a + ,b + )-S Q (a-,b-))P t , (27) 

where P t , It, Jt ar e taken from flip]). Further, it is easy to check that the 
following relations are true 

((a + ) T (b + )i - (6 + ) T (a+)i) - ((o_) T (6_), - (&-)r(a-)i) = 
d 



dl 
d_ 
dr 



(6+ 



6_) (fl+) - + (Q ' } - - (a + - a_) (M - + (Mt 



(a+h + (a_b 



(a+ — a_ 



(6+)j + (&_), 



p - ~rv ; 

ar 



u+ ((o+) T (6 + )i - (6+) T (a+)i) - u_ ((a_) T (6_)j - (6_) T (a_)j) 
((5 + ) T (6+), - (6+) T (5+)0 - ((5_) T (6_), - (6_) T (5_)0 = 



d 
~dl 
d 

dr 



(t+ _ t . ) (^). + (g-)-_ (s+ _ s . ) (M, + (t-), 



i - — « ; 
ar 



- K ((o+) T (6+)i " (6+)r(a+)i) " «- ((a-) T (6-)j - (&-)r(a_),)] 
((5+VCO+), - (a+) r (5' + )0 - ((£_) T (a_), - (a_) T (fiL.),) = 



d 
Jl 
d 



/ \ {S+)t + ('S'-)r / C 

(a + — a_) (o_| 



a- t 



(a+ — a_) 



+ esq, 



(5+ - s_ 



(a + )i + (a_); 



where Sk = S (ak,bk), (k = +,—)■ Taking into account these formulas 
relation fl27p produce ( ^5]) and inversely. 

□ 



REMARK 3.1 In case the time interval t is small and initial data are suffi- 
ciently smooth all terms in relation $2$ ) approximately equal to zero and one 
recovers formulas suggested in J7<|/ . Thus when the time elapsing from the 
moment of the emerging of a shock front is small the variational representa- 
tion is approximately valid. 
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REMARK 3.2 The nontrivial 2-D example which is known to the author 
where the variational representation works rigorously is the 2-D Riemann 
problem for (W with compressive piecewise constant potential velocity vec- 
tor. The example of corresponding initial potential So(a,b) is shown below 
(A,B,C are constants) 



S (a,b) = Ca- [B+ A ' ( 



j b ; b > max ^0, v^3aj 



So(a, b) = —Aa — Bb ; a > 0, V3a > b > —VSa 
So(a, b) = Ca + ^ 



B^jb ; b < min ^0, -VZa) , 



V3 

C > 0, A + V3B > 0, A - V3B > . 

This result is obtained when the discontinuity lines are straight and between 
every two neighboring lines is the same angle: 90° or 120° or 180°. The 
conjecture is that the angle value can be taken arbitrary. 

But there the another type of singularity arises: S-function in one point 
for density. 

The validity of variational principle in this case could be understood in 
the following way: the geometry of corresponding initial potential was flat 
and had high degree of symmetry. 

4 The flow description inside the shock. 

THEOREM 4.1 Suppose there exists the generalized solution to the problem 
(HP, /[TT\) in the form $12]), ftT3{), flTQ ) and the surface S can be defined from 
the parametric equations x = x(t,l), y = y(t,l). Then the following system 
of equations is true 

+ xi {V(g + - g_) - {q + v+ - Q-V-)} + 
yi{(Q+u+ - Q-uJ) -U(q+- Q-)} = 



tj + Xi {V(q + U + - Q-U-) - (q + U + V + - Q-U-V-)} + 



Vi { (Q+u 2 + ~ Q- u -) ~ u {Q+U+ - Q-U-)} = (21 



Jtj + xi {V(g + v + - Q-V-) - (g+vl - Q-V 2 _)} + 
Vi {i.Q+ u + v + ~ Q-U-V-) -U(q + v + - Q-vJ)} = 
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y = V , X = U , 

where U — It/Pt, V = Jt/Pt an d 'dot' denotes the differentiation with respect 
to t. 

PROOF. Taking into account (f|), (]5|) one has 



Further in this proof for our convenience we will do all calculations inde- 
pendently for indices "+", "— " and so omit indices in the expressions and 
write a, b, Qq, u, v, D instead of a i; b{, g l , u l , v$, D u i = +, -. From © 
one can infer after differentiation the relations for xi, yi, x, y. 

Now consider the right hand side of the first equation from (p8|) and write 
all terms concerning the index "+" (for "— " the calculations are analogous). 
One has 



x(t,l) = di + tu l (ai,bi) 
y(t,l) = bi + tv l (a, h bi) 
Qi = Qo(di,bi)/D(ai,bi) , 



(29) 



where (i = +,—) and 



D{a h bi) = + * K 



)J (l + *(«5 



g [xi (V -v) + y t (u -U)]= g [x t (y - v) + y\ (it - x)\ 
g [ai (1 + tu a ) + ktub] b (1 + tv b ) + atv a - 



g [aitv a + bi (1 + tvb)} a (1 + tu a ) + btiib 




Thus the first equation from ( p8|) is equivalent to the first equation from (p!6|). 
Taking into account the relations 




and 




one obtains the other two equations 



□ 



COROLLARY 4.1 The system @) is well defined. 
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PROOF. We have to prove that P t > under the dynamics of (28). Let 



us integrate P t with respect to / from some Iq to Iq + Al, where Al is small 
enough. Note that P t dl is exactly the sum of the areas which are 

bounded by the curves (i = +,—)'■ 

(oi(r, l ), bi(r, l )) , (^(r, l + Al), b^r, l + Al)) , < r < t ; 

(a,(t, 0, 0) , (Oi(0, 0> 0) > l < I < k + Al . 
Since A/ is arbitrary then P t > as t > 0. 

□ 

REMARK 4.1 The system ( p^ ) nonstrictly hyperbolic system with one 
eigenvalue and three eigenvectors. From the first three equations of system 

( ffify it is obvious to find the additional relation between (Pt^j , (it) > (jy 
by eliminating x\ and yi. 

REMARK 4.2 From the Cauchy-Kovalevskaya theorem one immediately gets 
the local existence and uniqueness theorem for (ffify in case of analytic coef- 
ficients and initial data. 

5 Specific type of motion. Constant external 
state. 

The system (E8J) does not satisfy Friedrichs' symmetrizability condition so 
to obtain the existence theorem encounters some problems. Nevertheless the 
internal dynamics of (|28|) is highly nontrivial which we demonstrate here in 
the simplest case of constant external density and velocity. 

Suppose that in ([28J) g + = const = g, g_ = const = g, u + = v + = 0, 
w_ = const = u, f_ = const = v; the velocity vector (u, v) satisfies condition 
II). To simplify the notations we also drop index 't' and 'waves' in (28|). Then 
one has 



(30) 



P+(g-g) {yiU - x t V} - g {yiu - x t v} = 
I + gu {yiU - xiV} - gu {y t u - x t v} = 
j + gv {ydJ - xiV} - gv {yiu - x t v} = 
x = U = I/P , y = V = J/P. 

From ( pCf ) one immediately gets the first integrals 

uJ - vl = u J - vl = C(l) . (31) 
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and 

uy-vx = ^-. (32) 

Now let us take the special initial conditions such that C(l) = 0. In other 
words there exists such k Q (l) that 

I = k u , J = k v . (33) 

Then taking into account fl3~l|) there exists some unknown k(t, I) such that 

I = ku , J = kv . 

And from (02D 



uyi - vxi = uy' - vx' = G{1) > . 

due to condition II). After rather simple transformations one arrives to the 
following system 

P + G(l) \k(g-g)-g) =0 

1 J (34) 



(Pk\ + gG(l) =0 



where k = k/P. 

Now expressing k from the first equation of and substituting it into 
the second equation one can find the expressions for P and k which read 

P 2 = p Q {lf - 2G(l)P (l)N(l)t + qqG{1)H 2 

P (35) 



k(g-g) = Q 



G(l) 



where N(l) = ko(l) (g — g) — g. 

From (|30| ) one obtains that x = ku, y = kv and the stability condition 
(see Definition 2) takes the form 



< J k{r,l)dr < t . (36) 
o 

It is easy to find that from (|36|) taking into account (j3~5D follows the 
stability condition (see Definition 2) 

< jfeo(0 < 1 • (37) 
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Then from one can easily see that N(l) < and so ( pop are well defined. 
From fl35p also follows that k — > x as t — > oo, where 

x^ ^ (38) 

Thus we have proved the following theorem 

THEOREM 5.1 Suppose that g + = const = g, g_ = const = g, u + = v + = 
0, u_ = const = u,v^ = const = v. 

Suppose also that G(l) > and ftSty is true. Than there exists the solution 
to the problem (fjp ; (fjjj and the shock front tends to the following one as 
t — > +oo 

x(l, t) = X (l) + KUt 

y(!,t) = yo(l) + mt , 

where x is taken from (PI 



Finally to illustrate the nontrivial character of the problem even with 
constant external fields let us derive the equation in the case C{1) ^ 0, but 
q = g. Then from ( |30l) one infers 

P = Q{yiu-xiv} 

+ \P~ U ) (39) 
±= P- 

Differentiating the first equation from (|3^|) with respect to t and taking into 



account integral (|32|) one gets the equation for P 

So even in the simplest case our system delivers us rather unusual equations 
of type (flop. To illustrate this let us perform the stability analysis for small 
perturbations to the model linear equation 

P = KP X , (41) 



where K = const. One have to find partial solutions to flil ) in the form 

Pit, x) = e ^ x+xt) , (42) 
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where i 2 = — 1, ( G R, A = c + iA, a G 1, A e 1. One immediately gets the 
restriction to the choice of £ and A 

A 2 = -Kii . 

Further, one has 

o 2 = A 2 ; 2aA = -K£ . (43) 
From (f43"D it follows that 

2A 2 = ±K£ , 

so we can choose such signs that —A ~ const^ and —A tends to +oo as 
£ tends to +oo. Thus in ( f4"2"D one has an arbitrary rapid growth of small 
perturbations with high frequencies and equation (|4l|) is ill-posed in the 
class of functions of finite smoothness because it does not satisfy classical 
Petrovsky condition. 



Appendix. 

Let us carry out some heuristic calculations to obtain the system (|28|) . The 
generalized solution ([]]), (fll"| ) can also be written in the form 

u = x, y) + (u + (t, x, y) - u_(t, x, y)) H(S) 

v = v _(t, x, y) + (v+(t, x, y) - «_(*, x, y)) K(S) (44) 

q = Q-(t, x, y) + (g + (t, x, y) - g_(t, x, y)) R(S) + A| 5=0 5(5) , 

where 5 is usual Dirac 5-function, but H, K, R are different Heaviside func- 
tions which can be distinguished by means of the following heuristic multi- 
plication formulas 

5 ■ H = st-8 , 5 ■ K = s 2 ■ 5 , (45) 

where si, s 2 are some functions on the surface S (in what follows we need not 
multiplication with R, so it is not included in (|4"5|)). In addition the following 
rather natural formulas are supposed to be true in the sense of distributions 

H 2 w K 2 ps HK w RH ps RK ta H . (46) 

Let us note that these formulas can be treated rigorously, for example, with 
the help of theory of new generalized functions 0, M (all basic ideas and 
lines in application to physics can also be found in 0) but here we do not 
need such rigor. 
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ties 



We can write 

U = U- + s\{u + — u ) , V = v _ + S2(v+ — f-) • 
Then taking into account the relations one has the following equali- 

QU = + (Q+U+ — Q-U-)H + XU5 

gv = + (q+v+ — Q-vJ)H + XV5 

gv 2 = g-u 2 _ + (g+u\ — g-u 2 _)H + XU 2 8 

guv = g^u^v_ + (g + u + v + — g_u_vJ)H + AC/V5 

^ 2 = g_ v 2 _ + _ q_ v 2 _)H + AV 2 <5 

Substituting these equalities in the system (|l|) and equating to zero expres- 
sions with different kind of singularities one obtains 

s t + us x + vs y = o 

\t + (XU) s + (XV) y +(Q+-Q-)S t + 
(g+u+ - g-u-)S x + {g+v + - g-V-)S y = 



(47) 



(XU) t + (XU 2 ) x + (XUV) y + {g + u+ - g.u-)St + 
(g+u 2 + — g-U 2 _)S x + (g + u + v + — Q-U-V-)S y = 

(XV) t + (XUV) X + (XV 2 ) y + (g + v + - g_v_)S t + 
(g + u + v + - g^u^v_)S x + (g+v+ - g-V 2 _)S y = . 

Now let us mention that the surface S can be defined from the equation 
S = x — X(t,y) = but the functions A, U, V depend only on t,y. Intro- 
ducing the differentiation along the direction (U, V) which will be denoted 
by 'dot' rewrite the system (f|7]) in the form 

y = V , X = U , y{0,l)=l 

X + XV y + X y (V(g + - g_) - (g+v + - Q-V-)) + 

(g+u + - q-U-) -U(g + - g_) = 



XU + X y [V (g+u + - g^u_ — U(g + — qJ)) + 
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U(g+V + - Q-V-) - (Q+U+V+ - Q-U-V-)] + 

U (U(g+ - qJ) - (g + u + - Q-uJ)) + 

{Q+u\ — Q-U 2 _) — U(q + U + — Q-U-) = 



(48) 



+ X y [V {q+v+ - q-V- - V{g + - + 
V(Q+v+ ~ Q-V-) - (q+v 2 + - Q~v 2 _)] + 
V(U(g + - g_) - (q+u+ - Q-U-)) + 
(g + u + v + — ^_?i_f_) — U(g + v + — Q-vJ) = . 

Now let us transform the system (f!8|) in the following way. Multiply 
the first equation by U and add to the second equation, then multiply the 
first equation by V and add to the third equation. Finally multiplying all 
obtained equations by yi and using the relation y\ = V y yi one gets the system 
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